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Abstract 

A^-Higgs doublet models (NHDM) are a popular framework to construct electroweak 
symmetry breaking mechanisms beyond the Standard model. Usually, one builds an 
NHDM scalar sector which is invariant under a certain symmetry group. Although 
several such groups have been used, no general analysis of symmetries possible in the 
NHDM scalar sector exists. Here, we make the first step towards this goal by classifying 
the elementary building blocks, namely the abelian symmetry groups, with a special 
emphasis on finite groups. We describe a strategy that identifies all abelian groups 
which are realizable as symmetry groups of the NHDM Higgs potential. We consider 
both the groups of Higgs-family transformations only and the groups which also contain 
generalized CP transformations. We illustrate this strategy with the examples of 3HDM 
and 4HDM and prove several statements for arbitrary N. 

1 Introduction 

The Higgs mechanism of the electroweak symmetry breaking (EWSB) has not yet been discov- 
ered experimentally, although the LHC is expected to provide soon first definite clues about 
the nature of the EWSB. Many different variants of the Higgs mechanism have been proposed 
so far [1]. Moreover, even if one focuses on a specific Higgs mechanism beyond the Standard 
Model (SM), one has many free parameters describing the scalar and Yukawa sectors, all of 
them having phenomenological consequences. Ideally, one would want to have a formalism 
that reconstructs the phenomenology of a given model in its most general form, that is, in 
the entire space of the allowed free parameters. Unfortunately, such a construction seems 
impossible for most non-minimal Higgs mechanisms. In this situation the standard approach 
is to narrow down the freedom by imposing certain symmetries both on the Higgs sector 
and on the Yukawa interactions. The problem is that in most cases it is not known a priori 
which symmetries can be used. Many specific symmetry groups have been analyzed but their 
complete lists for every class of non-minimal Higgs model are still lacking. 

The general task of classifying symmetries in models beyond SM can be, roughly, separated 
into two subtasks. First, one needs to know which symmetries can be implemented in the scalar 
sector of the model and how these symmetries are broken under EWSB. Second, one has to 
understand how these symmetry patterns can be encoded in the Yukawa sector of the model 
and what are the resulting properties of the fermions. 
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In this paper we exclusively focus on the first task in multi-Higgs-doublet models, which 
is an ambitious program on its own. It is true that the symmetries to be obtained within the 
scalar sector can be broken by a sufficiently complicated Yukawa sector. However, it does not 
diminish the importance of rigorous study within the scalar sector due to following reasons. 
(1) At the tree level, the vacuum properties in the electroweak theory come entirely from the 
scalar potential, |2]. Fermions can influence the vacuum properties only via loop corrections. 
Therefore, it is logical first to derive the tree-level features of the model, and next to see 
the effect of corrections. (2) The Yukawa couplings do not necessarily have to break every 
symmetry found in the scalar sector. It is easy to construct a Yukawa sector, for example by 
coupling all fermions to only one Higgs doublet, which inherits many of the symmetries of the 
scalar sector. This construction can lead to interesting phenomenological consequences, for 
example to the possibility of Zp-stabilized scalar dark matter, [3]. Which symmetry groups 
can be extended to the entire lagrangian is a separate question to study. (3) Historically, 
one of the most studied model beyond SM is the two-Higgs-doublet model (2HDM), [11[5]. 
After decades of theoretical and phenomenological studies which included various forms of 
the Yukawa sector, it became clear that there were several pressing questions (including the 
origin of CP-violation in 2HDM) which required an in-depth look into the properties of the 
general scalar potential. Only when sufficient progress in the scalar sector was made, did it 
become clear how to construct Yukawa sectors incorporating additional CP-violating features, 
see examples in [6]. 

1.1 Realizable symmetries 

When searching for the symmetry group which can be implemented in the scalar sector of a 
non-minimal Higgs model, one must distinguish between realizable and n on- realizable groups. 
If it is possible to write a potential which is symmetric under a group G but not symmetric un- 
der any larger symmetry group containing it, we call G a realizable group (the exact definition 
will be given in Section [273|1 . A non-realizable group, on the contrary, automatically leads to 
a larger symmetry group of the potential. This means that if we write a potential symmetric 
under a non-realizable G, we'll discover that it is in fact symmetric under a larger symmetry 
group if D G, and there is no way to avoid that. Roughly speaking, it means that the Higgs 
potential of a given model cannot accommodate overly complicated symmetry patterns. 

Thus, the true symmetry properties of the potentials are reflected in realizable groups. 

Just to give an example, cyclic groups such as Zy and Z13 have been used in the context 
of the flavor symmetry problem, see e.g. |^ and references therein. These symmetry groups 
are usually implemented in the space of fermion families, but sometimes they are transferred 
to the Higgs sector of the model. Anticipating the results of this paper, we can state that 
these symmetry groups are not realizable in the three-Higgs-doublet model potential: trying 
to impose such a symmetry, one ends up with a potential that has a continuous symmetry. 

Therefore, if one focuses on the scalar sector of the model and aims at classifying the 
possible groups which can be implemented in the scalar potential, it is natural to restrict 
one's choice to the realizable symmetry groups. 
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1.2 Multi-Higgs-doublet models 

One particular non-minimal Higgs sector where the classification of possible symmetries of 
the scalar sector has been established is the two-Higgs-doublet (2HDM), [HHUH]. Although 
the general Higgs potential in 2HDM cannot be minimized explicitly with the straightforward 
algebra, it can still be analyzed fairly completely with several basis-invariant tools that were 
developed in the last years. The basic idea is that reparametrization transformations, or the 
Higgs-basis changes, link differently looking but physically equivalent potentials to each other. 
This idea can be implemented via the tensorial formalism at the level of Higgs fields, [HHTT] , or 
via geometric constructions in the space of gauge-invariant bilinears, [T214T7] . In particular, all 
possible groups of Higgs-family symmetries and generalized CP-transformations of the 2HDM 
potential have been found. In addition, it was noted that if one disregards the U{1) gauge 
transformations the structure of the Higgs potential of 2HDM is invariant under a wider class 
of transformations, |18] . 

It is natural to attempt to extend these successful techniques to N Higgs doublets. Some 
properties of the general NHDM potential were analyzed in [T ^[TU] - [2^ . In particular, the 
bilinear-space geometric approach was recently adapted to the NHDM in [25l[26]. In what 
concerns the symmetries of the scalar sector of NHDM, several groups have attacked various 
issues of this problem, [22142^ I261429] . In particular, in [22] an accurate and lengthy analysis 
of possible textures of the quartic interaction in the 3HDM led to the list of one-parametric 
symmetries of the potential. Although the method of [22] could be extended to > 3, its 
practical implementation relies on the set of all possible textures for the quartic coupling 
constants, which is not known a priori. 

In [26] an attempt was made to describe possible symmetry groups of NHDM in the adjoint 
representation. The symmetry properties of the vectors Mj and Aqj and the symmetric tensor 
Ajj in the expression (jl]) written below could be easily identified, but the main difficulty was 
to take into account the non-trivial shape of the NHDM orbit space, see details in Section 
12.11 Mathematically, this led to the problem of classifying intersections of two orthogonal 
subgroups of SO{N^ ~ 1)) which was not solved in [26] . 

It is therefore fair to state that despite several partial successes the task of classifying 
all possible Higgs-family and generalized- CP symmetries in NHDM is far from being accom- 
plished. 

In this paper we make a step towards solving this problem. We focus on abelian symmetries 
of the NHDM scalar potential and describe a strategy that algorithmically identifies realizable 
abelian symmetry groups for any given A^. The strategy first addresses groups of unitary 
transformations in the space of Higgs families and then it is extended to groups that include 
antiunitary (generalized CP) transformations. The strategy also yields explicit examples of 
the potentials symmetric under a given realizable group. 

The structure of the paper is the following. In Section |2] we describe the Higgs potential 
of the general NHDM, discuss the group of reparametrization transformations and introduce 
the notion of realizable symmetries. In Section [3] we describe the strategy that identifies all 
realizable abelian symmetry groups. Then in Section H] we illustrate the general approach 
with the examples of 3HDM and 4HDM. Section [S] contains some further general results valid 
for all A^. Finally, in Section [6] we make several remarks and draw our conclusions. In the 
Appendices we prove the theorem about maximal abelian groups of PSU{N) and work out in 
detail the groups containing antiunitary transformations in the case of 3HDM. 
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2 Scalar sector of the A^-Higgs doublet model 



2.1 NHDM potential 

In the NHDM we introduce complex Higgs doublets with electroweak isospin Y = 1/2: 



a 
4J) 



l,...,iV. (1) 



Hermitian conjugation is denoted by dagger: 0]^. The generic renormalizable Higgs potential 
in NHDM can be written in the tensorial form [Il[2]: 

V = Ya,{<P%) + ZaU<P%M<Pd) , (2) 

where all indices run from 1 to A^. Coefficients in the quadratic and quartic parts of the 
potential are grouped into components of tensors Yab and Zabcd, respectively; there are A^^ 
independent components in Y and N^lN"^ + l)/2 independent components in Z. 

Once the Higgs potential is given, the first task is to find its global minimum. To this 
end, we replace the Higgs field operators by their vacuum expectation values, and interprete 
([2]) as a scalar function defined in the 2A^-complex-dimensional vector space C^^. The elec- 
troweak gauge group is reduced to the global SU{2) x f/(l) transformation group acting in 
this space. Pick up a point x G C^^ and apply a transformation g G SU{2) x U{1). All 
points = g ■ X which can be reached in this way form a gauge orbit. Since the potential is 
electroweak-symmetric by construction, its values are equal at all points throughout a chosen 
orbit. Therefore, the Higgs potential can be throught of as a scalar function defined in this 
orbit space. 

This leads to the question of how this orbit space can be described. In principle, one can 
uniquely define an orbit by providing A^^ gauge-invariant bilinears not all of them 

being algebraically independent. An even more convenient way is to introduce the real vector 
= (ro, Ti) with 

= \/ ~2N~ '^'^'^^ ' ^ \^ii^i)ab^b , (3) 

where Aj are the standard SU{N) generators. Then, the orbit space can be described com- 
pletely by a list of A^ algebraic equalities and inequalities written for tq and rj with coefficients 
being the invariant SU{N) tensors of rank < A^, [25]. The Higgs potential can be rewritten 
as ^ ^ 

V = -Mora - Min + -Aooro + Aoiror^ + -Kinrj , (4) 

which, if needed, can be represented even more compactly a.sV = —M^r^ + K^^r^r'^ /2. All 
the coefficients Yah and Zabcd become now components of the vector and symmetric tensor 
A^jy, which transform under adjoint representation of the group GL{N,C). For more details 
about the properties of the Higgs potential in the space of bilinears, see [26] . 



2.2 The group of reparametrization transformations 

When discussing symmetries of the potential, we focus on the reparametrization transfor- 
mations, which are non-degenerate linear transformations mixing different doublets (pa but 
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keeping invariant the kinetic term (which includes interaction of the Higgs fields with the 
gauge sector of the model). Alternatively, they can be defined as norm-preserving transfor- 
mations of doublets that do not change the intradoublet structure. In this work we do not 
use the more general transformations in the spirit of [H] but focus on the "classically defined" 
reparametrization transformations . 

A reparametrization transformation must be unitary (a Higgs-family transformation) or 
antiunitary (a generalized CP-transformation): 

U : (pa^ Uab(t)h or Ucp = U ■ J : 0^ ^ ^afo^I , (5) 

with a unitary matrix Uab- The transformation J = CP acts on doublets by complex conju- 
gation and satisfies = 1. 

Let us focus first on the unitary transformations U. A priori, such transformations form the 
group U{N). However, the overall phase factor multiplication is already taken into account by 
the U{1)y from the gauge group. This leaves us with the SU{N) group of reparametrization 
transformations. Then, this group has a non-trivial center Z{SU{N)) = generated by the 
diagonal matrix exp(27rz/A^) ■ In, where Iat is the identity matrix. Therefore, the group of 
physically distinct unitary reparametrization transformations is the projective special unitary 
group 

Gu = PSU{N) ~ SU{N)/Zn . (6) 

Consider now antiunitary transformations Ucp = U ■ J , U E U{N). One can define an 
action of J on the group U {N) by 

J -U -J = U\ (7) 

where asterisk denoted complex conjugation. This action leaves invariant both the overall 
phase subgroup U{1)y and the center of the SU{N). Therefore, we again can consider only 
Ucp = U ■ J such that U G PSU (N). Once this action is defined, we can represent all distinct 
reparametrization transformations as a semidirect product 

G = PSU{N) y<i Z2 . (8) 

Recall that G is a semidirect product of subgroups A and B with normal subgroup A (this 
fact is usually denoted as G = A xi B) if A is normal in G and the intersection AilB is trivial 
(equals {1}). 

2.3 Realizable symmetry groups 

A (reparametrization) symmetry of the Higgs lagrangian is such a reparametrization trans- 
formation that leaves the potential invariant: V{U(^cp)4') = ^(0)- For any given potential 
all such transformations form a group, the automorphism group of the potential, which we 
denote as Gy = Aut(y) fl G and which is obviously a subgroup of G in (jHl). 

Let us stress that when we say a group Gy is a symmetry group of the potential V, we 
mean that it is equal to, not just contained in, Aut{V) fl G. In [26] it was suggested to call 
such groups realizable symmetry groups. In simple words, a symmetry group Gy is realizable 
if one can write a potential symmetric under Gy but not symmetric under any overgroup of 
Gy, i.e., over H with Gy S ^ < G. For example, a potential which depends on the first 
doublet only via (0l0i) is obviously symmetric under the cyclic group Z„ of discrete phase 
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rotations of this doublet for any n. However, these Z„'s have no interest on their own because 
they trivially arise as subgroups of the larger symmetry group of this potential U{1) describing 
arbitrary phase rotations. It is this U{1), not its individual subgroups Z„, which has a chance 
to be the realizable symmetry group of this potential. 

Since in this paper we deal with groups Au containing only unitary transformations of 
Higgs families and with groups A that can contain antiunitary transformations, let us adopt 
the following convention: a group A^ G is called realizable if A^ = Aut(y) fl and a 
group A G G is called realizable ii A = Aut{V) fl G. 

Note that in the case of 2HDM, where all possible reparametrization symmetries are known, 
Z2, (2,2)^, (2,2)^, f/(l), U{1) X Z2, SU{2), each of these symmetry groups is realizable, which 
is most easily proved in the space of bilinears. 

3 Finding abelian groups 

We now focus on abelian groups for NHDM and describe in this Section the strategy that 
identifies all realizable abelian groups of the NHDM potential for any A^. This strategy can be 
outlined as follows: we first describe maximal abelian subgroups of PSU{N), then we explore 
their realizable subgroups, and finally we check which of these groups can be extended by 
including antiunitary transformations. 

3.1 Heuristic explanations 

Before presenting the rigorous strategy which characterizes realizable abelian groups in NHDM, 
we first give heuristic arguments which should facilitate understanding of the main idea. 

Suppose we are given a potential V of the A^-doublet model and we want to find which 
phase rotations leave this potential invariant. Clearly, these phase rotations form a group 
which is a subgroup of the group of all possible phase rotations of doublets, that is, the 
group of all diagonal unitary N x N matrices acting in the space of doublets. This group is 
[[/(l)]^ C U{N) and can be parametrized by parameters aj G [0,27r): 

diag[e*"\ e*'^^ . . . , e^"^] . (9) 

The potential is a collection of k monomial terms each of the form {(^^(pb) or {4'\4'b){4'l4'd)- 
Upon a generic phase rotation dH]), each monomial term gains its own phase rotation. For 
example, (0^06) with a b gains the phase at — aa, (0i06)(0l0c) with a, b, c all distinct gains 
the phase + — '^c^a, etc. In short, each monomial gets a phase rotation which is a linear 
function of a's with integer coefficients: Ylf=i^j^j- The vector of coefficients mj can be 
brought by permutation and overall sign change to one of the following forms: 

(1,-1,0,...), (2,-2,0,...), (2,-1,-1,...), (1,1,-1,-1,...), (10) 

or a zero vector. Note that in all cases Ylj '^j = 0- Thus, the phase transformation properties 
of a given monomial are fully described by its vector rrij. The phase transformation prop- 
erties of the potential V , which is a coUeciton of k monomials, is characterized by k vectors 
"^ij, ''^2j, • • • , ■'^fcj, each rrtij being of one of the types in ( fTOj) . 

If we want a monomial to be invariant under a given transformation defined by phases {oj}, 
we require that Ylj=i''^j'^j — with some integer n. If we want the entire potential to be 
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invariant under a given phase transformation, we require this for each individual monomial. 
In other words, we require that there exist k integers {rij} such that the phases {aj} satisfy 
the following system of linear equations: 

N 

^^mij-aj = 27ini , for all 1 < « < /c . (11) 

i=i 

Solving this system for {aj} yields the phase rotations that leave the given potential V in- 
variant. 

One class of solutions can be easily identified: if all aj are equal, aj = a, then ( ITTll with 
TT-i = is satisfied for any a. These solutions form the U{1) subgroup inside [[/(l)]^ and simply 
reflect the fact that the potential is constructed from bilinears </>;,). These solutions become 
trivial when we pass from the U (N) to the SU{N) group of transformations. However, there 
can exist additional solutions of f|TT]) . They form a group which remains non-trivial once we 
pass from U{N) to SU{N) and further to PSU{N), which is the group of physically distinct 
Higgs-family reparametrization transformations. It is these solutions which we are interested 
in. 

In order to find these solutions, we note that a matrix with integer entries can be "diag- 
onalized" by a sequence of elementary operations on its rows or columns: permutation, sign 
change, and addition of one row or column to another row (column). "Diagonalization" for 
a non-square matrix means that the only entries rriij that can remain no n- zero are at z = j. 
After that, the system splits into k equations on phases of the generic form 

rui^iai = 27ihi , tt^ G [0, 27r) , Ui e Z , (12) 

with non- negative integer rrii^i. If rrii^i = 0, this equation has solution for any af, thus i-th 
equation gives a factor Z to the symmetry group of the potential. If rrii^i = 1, then this equation 
has no non-trivial solution, and the i-th equation does not contribute to the symmetry group. 
If rrii^i = di > 1, then this equation has di solutions which are multiples of at = 27r /di, and 
the i-th equation contributes the factor Z^- to the symmetry group. The full symmetry group 
of phase rotations is then constructed as direct product of these factors. 

Thus, the task reduces to studying which diagonal values of the matrix rrii j can arise in 
a model with doublets. For small values of A^ this task can be solved explicitly, while for 
general A^ one must rely upon subtle properties of rriij which stem from ffTOl) . 

3.2 Maximal abelian subgroups of PSU{N) 

In the previous subsection we outlined the main idea of the classification strategy. However 
we worked there in the group U{N), while the reparametrization group is PSU{N), which 
essentially complicated the analysis. Repeating this strategy for abelian subgroups of PSU (N) 
is the subject of the remaining part of the text. 

We start by reminding the reader of the definition of a maximal abelian subgroup. A 
maximal abelian subgroup of Gu ([6]) is an abelian group that is not contained in a larger 
abelian subgroup of Gu- A priori, there can be several maximal abelian subgroups in a given 
group. Any subgroup of Gu must be either a maximal one, or lie inside a maximal one. 
Therefore, we first need to identify all maximal abelian subgroups of Gu and then study their 
realizable subgroups. 



7 



If Gu were SU{N), then the situation would be simple. As we describe in Appendix A, all 
maximal abelian subgroups of SU (N) are the so called maximal tori, 

[f/(l)]^-i = f/(l) X f/(l) X ■ ■ ■ X f/(l) . (13) 

All maximal tori are conjugate inside SU{N), that is, given two maximal tori Ti and T2 there 
exists g e SU{3) such that g~^Tig = T2. It means that without loss of generality one could 
pick up a specific maximal torus, for example, the one that is represented by phase rotations 
of individual doublets 

diag[e*'^S e*"^ . . . , e^"^-\ e"^^"'] , (14) 

and study its subgroups. The analysis would then proceed essentially as we explained in the 
previous subsection with an additional condition that all a's sum to zero. 

However, the group of distinct reparametrizations is Gu = PSU{N), and it has a richer 
structure. Consider the canonical epimorphism (i.e., surjective homomorphism) 

~: SU{N) SU{N)/Z{SU{N)) = SU{N) ~ PSU{N) , (15) 

where Z{SU{N)) ~ Z^v is the center of SU{N). Given subgroup A of PSU{N) we denote by 
A G SU{N) the complete preimage of A, i.e., A = {g E SU{N) \ g E A}. If A is abelian, its 
complete preimage A does not have to be abelian. Specifically, in Appendix [A] we prove the 
following Theorem: 

Theorem 1. Let A be a maximal abelian subgroup of PSU{N), denote by A the complete 
preimage of A in SU (N) . Then one of the following holds: 

(1) A is abelian and A is conjugate to the subgroup of all diagonal matrices in SU{N). 

(2) A is a finite nilpotent group of class 2 (i.e., the commutator [A, A] lies in the cen- 
ter Z{A)), the exponent of A is divisible by N and divides N"^ . Moreover, Z{A) = 
Z{SU{N)). 

This theorem states that there are two sorts of maximal abelian groups in PSU{N): (1) 
maximal tori, which will be constructed below, (2) certain finite abelian groups, which are not 
subgroups of maximal tori and which must be treated separately. 

3.3 Maximal tori in PSU{N) 

Let us explicitly construct the maximal torus in PSU{N) referred to in Theorem 1. We first 
introduce some convenient notation. A diagonal unitary matrix acting in the space of Higgs 
doublets and performing phase rotations of individual doublets, such as f|T^ . will be written 
as a vector of phases: 

(ai, a2, . . . , aN-i, - • (16) 

In addition, if Mi, . . . , are subsets of a group G, then (Mi, . . . , M^) denotes the subgroup 
generated by Mi U . . . U M^ of G. Then we construct a maximal torus in SU{N) 

To = (f/(l)i,t/(l)2,--- ,t/(l)^-i), (17) 
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where 



f/(l)i = 1, 0, 0, 

t/(l)2 = «2(-2, 1, 1, 0, 
t/(l)3 = a3(-3, 1, 1, 1, 

f/(l)^_i = a^_i(-iV+l, 1, 1, 1, 1), (18) 

with all ai G [0, 27r). Clearly, for every j we have (f/(l)i, . . . , [/(l)j_i) fl f/(l)j = {e}, whence 
To is a direct product of f/(l)i, . . . , f/(l)7v-i- In particular, any element u G SU{N) can be 
uniquely written as 

M = Ui(ai)u2(a2) ■ ■ ■MAr-i(a7V-i) , UieU{l)i. (19) 

Moreover, the center Z{SU{N)) is contained in the last group and is generated by ctAr-i = 
2ti/N. One can therefore introduce t/(l)^_^ = t/(l)Ar_i/Z(S'?7(A^)), which can be parametrized 
as 

—r^ / iV- 1 1 1 \ 

where a^-i G [0, 27r). Since Z{SU{N)) n f/(l)j = 1 for j = 1, . . . , - 2, we identify these 
subgroups with their homomorphic images and omit the bar. Thus the homomorphic image 
T of To can be written as 

T = f/(l)ixt/(l)2X...xZ7(iy^_,. (21) 
3.4 Identifying symmetries of the potential 

Next we study which subgroups of the maximal torus T can be realizable in the scalar sector 
of NHDM. 

We start from the most general T-symmetric potential: 

V{T) = -Y.ml{<Pi<Pa) + Y^X,,[<pl<p,)[<pl<p,) + Y,Kb{<P%M<t^a) , (22) 

a a,b aj^b 

where among the N{N — l)/2 terms in the last sum there are only 2N — 3 algebraically 
independent ones, [25]. Each term in this potential transforms trivially under the entire T. 
The important fact now is that a sufficiently general potential of this form has no other unitary 
symmetry. In fact, an even stronger statement is true: 

Theorem 2. Consider potential V{T), l[2^) . defined by the set of coefficients m\, Xab, A'^j,. 

1. There exist coefficients such that the only unitary symmetries of the potential are the 
phase rotations from T. 

2. If one constructs V = V(T) + Vi by adding further terms that were absent in ^ME}> then 
for any Vi there still exist coefficients of V{T) such that any unitary symmetry of V 
belongs to T . 



...,0), 
...,0), 
...,0), 
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Proof. The first assertion is a particular case of the second, which we now prove. Consider 
a generic potential V in the space of bilinears (jlj). Any unitary transformation of the Higgs 
families corresponds to an orthogonal transformation in the rj-space. Therefore, any uni- 
tary symmetry of the potential must leave invariant each of the five terms in (jlj). Consider, 
specifically, the second and the forth terms: 

V = -Ya,{<P%) + ZaM<Pb) 5^(0l0c) + . . . , (23) 

c 

where we omitted the remaining terms and introduced the tensor 

Zab = Zabdd — J^^ab Z^cdd i (24) 

d c,d 

which is basically the of |9] made traceless. Both Yat and Zab are hermitean matrices, and 
their diagonal values can be arbitrarily changed by changing and Xab in fl22p . In particular, 
one can adjust them so that the spectra of Yat and Zab are non-degenerate. Their eigenvectors 
can be then used to construct two orthonormal bases in the space of Higgs doublets, {e^-^''} and 
{e^-^''}, respectively {i = 1, . . . , N). Let us also introduce the canonical basis {e^^"'*} consisting 
of vectors of the form (0, . . . , 1, . . . , 0), so that the maximal torus T consists of phase rotations 
of these vectors. 

The symmetry group of the first term in (125]) is the maximal torus Ty of phase rotations of 
individual eigenvectors ef^\ The symmetry group of the second term in is the maximal 
torus Tz of phase rotations of individual eigenvectors e\ . The symmetry group of fl23|) is. 



therefore, Ty fl Tz, the phase rotations of common eigenvectors of the two matrices (and the 
symmetry group of the entire V is its subgroup). 

Take an eigenvector of Yah- Its orientation can be changed by readjusting coefficients if 
and only if it does not belong to the canonical basis. The same applies to Zab- Therefore, even 
if Yab and Zab had common eigenvectors not belonging to the canonical basis, one could avoid 
this coincidence by readjusting the coefficients and Xab- After this procedure, the only 
common eigenvectors of Yab and Zab will be those belonging to the canonical basis. Therefore, 
the symmetry group of ( 12^ . and consequently, of the full potential will be a subgroup of 
T. □ 

The theorem just proved guarantees that when we start from the T-symmetric potential 
(122 p and add more terms, we will never generate any new unitary symmetry that was not 
already present in T. This is the crucial step in proving that the groups described below are 
realizable. Our task now is to find which subgroups of T can be obtained in this way. 

Consider a bilinear (^^(pb, a ^ b. It gains a phase change under T (|2T]) which linearly 
depends on the angles af. 

(j)l(j)b -> exp[i{pabai + qaba2 H h tabaN-i)] ■ (t>i(t>b (25) 

with some integer coefficients pab, <lab, - - - , tab- Note that all coefficients are antisymmetric in 
their indices: pba = —Pab, etc. These coefficients can be represented by real antisymmetric 
matrices with integer values, or graphically, as labels of the edges of — 1 oriented graphs, 
one for each U{1) group. Each such graph has A^ vertices, corresponding to doublets (pa', all 
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vertices are joined with arbitrarily oriented edged, orientation indicated by an arrow. An edge 
oriented from (ph towards 0^ (edge 6 — )■ a) is associated with the bihnear (f>l(f>b and is labeled 
by pba in the first graph, qi,a in the second graph, etc. Examples of these graphs will be shown 
later for 3HDM and 4HDM. 

The Higgs potential is a sum of monomial terms which are linear or quadratic in 4>l^4>b- 
Consider one such term and calculate its coefficients p, . . . ,t. Let us first focus on how this 
term depends on any single U{l)i subgroup of T. There are two possibilities depending on 
the value of the i-th coefficient: 

• If the coefficient k in front of is zero, this terms is f/(l)j-symmetric. 

• If the coefficient k ^ 0, then this term is symmetric under the Z\k\ subgroup of U{l)i- 
group generated by phase rotations by 2TT/\k\. 

However, even if a given monomial happens to have finite symmetry groups with respect to 
each single f/(l)j, its symmetry group under the entire T is still continuous: for any set of 
coefficients one can adjust angles in such a way that Patc^i + ■ ■ ■ + tabdN-i = 0. Therefore, 
when studying symmetries of a given term or a sum of terms, we cannot limit ourselves to 
individual U{l)i groups but must consider the full maximal torus. 

The strategy presented below guarantees that we find all possible realizable subgroups of 
the maximal torus T, both finite and infinite. 

Consider a Higgs potential V which, in addition to the T-symmetric part (I22p contains 
k > 1 additional terms, with coefficients pi, qi, . . .ti to p^, g^, . . .tk- This potential defines 
the following (A^ — 1) x /c matrix of coefficients: 



X{V) 



( Pi qi 
P2 q2 

\ Pk qk 



tk J 



( 



\ rrik^i 



(26) 



mk,N-i J 



Here the second form of the matrix agrees with the notation of (ITTl) . The symmetry group 
of this potential can be derived from the set of non-trivial solutions for of the following 
equations: 

/ «i \ / 27mi \ 



XiV) 



\ aN-i J 



(27) 



There are two major possibilities depending on the rank of this matrix. 



Finite symmetry group. If rajakXiV) = N — 1, then there is no non-trivial solution 
of the equation fl27j) with the trivial right-hand side (i.e. all Ui = 0). Instead, there 
exists a unique solution for any non-trivial set of rij, and all such solutions form the 
finite group of phase rotations of the given potential. 

To find which symmetry groups can be obtained in this way, we take exactly — 1 
monomials, so that the matrix X{V) becomes a square matrix with a non-zero deter- 
minant. It is known that any square matrix with integer entries can be diagonalized 
by a sequence of the following operations: swapping two rows or two columns, adding a 
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row (a column) to another row (column) and multiplying a row (a column) by —1. Af- 
ter diagonalization, the matrix X(y) becomes diag((ii, . . . , d]\f^i), where di are positive 
non-zero integers. This matrix still defines the equation f l27|) for a^, which are linear 
combinations of a^'s and with n[ G Z. Therefore, the finite symmetry group of this 
matrix is Zdi x ■ ■ ■ x '^d^^i (where Zi means no non-trivial symmetry). 

Note also that each of the allowed manipulations conserves the absolute value of the 
determinant of X{V). Therefore, even before diagonalization one can calculate the 
order of the finite symmetry group as | detX(\^)|. 

This derivation leads us to the strategy that identifies all finite subgroups of torus re- 
alizable as symmetry groups of the Higgs potential in NHDM: write down all possible 
monomials with A^-doublets, consider all possible subsets with exactly — 1 distinct 
monomials, construct the matrix X for this subset and find its symmetry group follow- 
ing the above scheme. Although this strategy is far from being optimal, its algorithmic 
nature allows it to be easily implemented in a machine code. 

• Continuous symmetry group. If rankX(y) < A^— 1, so that D = A^— 1— rankX(V") > 
0, then there exists a dimensional subspace in the space of angles a^, which solves 
the equation fl27|) for the trivial right-hand side. One can then focus on the orthogonal 
complement of this subspace, where no non-trivial solution of the homogeneous equation 
is possible, and repeat the above strategy to find the finite symmetry group Gd in this 
subspace. The symmetry group of the potential is then [f/(l)]^ x Go- 

3.5 Identifying symmetries using the dual group 

The strategy just exposed can be put on a firmer algebraic ground with the classic formalism 
of the Pontryagin duality, [30] . 

Let T be a subgroup of all diagonal matrices of PSU{N), i.e., T is a maximal torus, 
and consider the decomposition T = Ui x . . . x Un, where n = N — 1. By construction, 
Ui = {ui{ai) I tti G [0,27r)}. By K we denote the factor group M/Z (we use additive notation 
for the operation in this group) and, given an abelian group A the homomorphism x : A — )■ A' 
is called a character of A. The set of all characters of A become an abelian group under the 
operation given by (xi + X2){o,) = Xi(o) + X2(a) and we denote this group by X{A). This is 
the dual group of A. If X{A) is found, then A can be recovered by considering X{X{A)) the 
dual group of X{A). 

Consider Xi ^ X(T) defined on generators by 



Clearly the cyclic group Xi := {xi) is isomorphic to X{Ui) and X{T) = Xi x . . . x Xn, i.e., 
Xi,---,Xn are free generators of a free abelian group X{T). Now consider a potential V, 
denote X{V) by M (recall that M is a /c x n- matrix with integer items) and denote the items 
of M by rriij. An element x = ui(ai) ■ . . . ■ M„(a„) G T lies in Aut(y) if and only if f l27|) 
holds. That is, x G Aut(V^) if and only if for each i we have C,i{x) = 0, where the character C,i 
is defined by = YJj=i ^hjXj- 
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Consider a subgroup H = (^i, . . . , ^fc) of X{T). The above arguments imply that H is the 
annulator of Aut{V), and that X{T)/H ~ X(Aut(\^)). If we diagonahze M by the standard 
transformations we obtain the following matrix 



D 



f di . . . \ 
4 ... 



V ... dk ... / 

where divides di for each z. Now [311 Theorem 8.1.1] implies that 

X{T)/H ~ Zrf, X . . . X Zrf, X Z"-^ 

where Z^. is a cyclic group of order di and Z"~'^ is a direct product oi n — k copies of Z. In 
view of the Pontryagin duality between groups and characters, we obtain that 

Aut(\/) ^Zrf, X ... X Zrf, X ?7, 

where U is an (n — /c)- dimensional torus. 



3.6 Antiunitary transformations 

Having found the list of realizable subgroups of torus T, we can ask whether these groups can 
be extended to larger abelian groups that would include not only unitary but also antiunitary 
transformations. Here we describe the strategy that solves this problem. 

As described above, one can define the action of the CP-transformation J on PSU{N) 
given by ([7]): J ■ U ■ J = U*. This action can be restricted to the maximal torus T ( |2T1) and 
even further, to any subgroup A C T (indeed, if a G A, then a* E A). However, J does not 
commute with a generic element of T, therefore the group {A, J) is not, in general, abelian. 

In order to embed A in a larger abelian group, we must find an antiunitary transformation 
J' = b ■ J which commutes with any element a E A: 

J'a{J')-^ = a ^ bJaJb'^ = a ^ ba~^b~^ = a ^ b = aba. (28) 

The last expression here is a linear matrix equation for the matrix b at any given a. Since a 
is diagonal, {aba)ij = aaajjbij, so whenever aaajj ^ 1, one must set bij to zero. 

If at least one matrix b satisfying (128|) is found, all the other matrices can be constructed 
with the help of the last form of this equation. Suppose b' = xb also satisfies this equation; 
then, axa~^ = x. Therefore, x can be any unitary matrix from PSU (N) commuting with a, 
that is, X centralizes the chosen abelian group A. 

In summary, the strategy for embedding a given abelian group A (ZT into a larger abelian 
group Acp which includes antiunitary transformations proceeds in five steps: 

• Find one b which solves the matrix equation (128 p for all a E A. 

• Find all x G PSU{N) which commute with each element of A, that is, construct the 
centralizer of A in PSU{N). 
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• Find restrictions placed on x's by the requirement that product of two different antiu- 
nitary transformations {xiJ') and {X2J') belongs to A. This procedure can result in 
several distinct groups Ac p. 

• Find which among the A-symmetric terms in the potential are also symmetric under 
some of Aqp] drop terms which are not. 

• Check that the resulting collection of terms is not symmetric under a larger group of 
unitary transformations, that is, A is still a realizable symmetry group of this collection 
of terms. 

As an exercise, let us apply this strategy to the full maximal torus T. A generic a E T 
acting on a doublet (pi generates a non-trivial phase rotation ipi{ai, . . . , aAr_i) which can be 
reconstructed from f|T8l) . The equation fl28|) then becomes 

eiii'^+i'j)bij = bij . (29) 

Since ipi + ipj ^ for any i,j, the only solution to this equation is bij = for all i and 
j. This means that there is no J' that would commute with every element of the torus T. 
Thus, T cannot be embedded into a larger abelian group Tcp that would include antiunitary 
transformations. 

In order to avoid possible confusion, we note that in fact J is a. symmetry of the T- 
symmetric potential. However, the symmetry group, generated by T and J, is a semidirect 
product T X J, and it is nonabehan. 

4 Examples of 3HDM and 4HDM 

In this Section we illustrate the general strategy presented above with the particular cases of 
3HDM and 4HDM. For 3HDM we give the full list of abelian groups realizable as symmetry 
groups of the Higgs potential and show explicit examples of such potentials; for 4HDM we do 
it only for the finite groups of unitary transformations. 

4.1 Maximal torus for 3HDM and its subgroups 

In the 3HDM the representative maximal torus T C PSU (3) is parametrized as 

T = t/(l)ixt/(l)2, t/(l)i = a(-l, 1, 0), t/(l)2 = /3 (^-| 1 1^ , (30) 

where a,f3 G [0,27r). The most general Higgs potential symmetric under T is 

V = -ml{4<P^) - ml{4<P2) - ml{4<f)s) (31) 

+ Ai2(0l0l)(0^02) + X23{4h){4h) + Ai3(0^03)(0l0l) 
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There are six bilinear combinations of doublets transforming non-trivially under T. The way 
they transform under f/(l)i and f/(l)2 is described by two coefficients p and q defined in fl25|) . 
Here is the list of these coefficients for the three independent bilinears: 





P 






-2 


-1 


(0^2) 


1 





(0l03) 


1 


1 



(32) 



For the conjugate bilinears, the coefficients p and q carry the opposite signs with respect to 
fl5^ . These coefficients are shown graphically as labels of the edges of two oriented graphs 
shown in Fig. [H 




Figure 1: Coefficients p and q as labels of triangles representing bilinears in 3HDM 

Before applying the general strategy described in the previous Section, let us check which 
finite symmetry groups can arise as subgroups of either f/(l)i or U{1)2 individually. 

Consider first f/(l)i. The coefficient p of any monomial can be obtained either by picking 
the labels from the first graph in Fig. [1] directly or by summing two such labels, multiplied by 
— 1 when needed. In this way we can obtain any \p\ from to 4. For a monomial with p = 0, 
the symmetry group is the entire U{l)i. For \p\ = 1, there is no non-trivial symmetry. For 
\p\ = 2,3,4, we obtain the cyclic group Z\p\. In each case it is straightforward to construct 
the monomials with a given symmetry; for example, (0i03)(0203) conjugate are U{l)i- 

symmetric, while (0^02) (^i^s), (0i02)(0302) and their conjugates are symmetric under the 
Zs-group with a generator 

a = {co, 1, uj^'^) , CO = exp (27rz/3) . (33) 

In the case of the group f/(l)2, the labels can sum up to |g| = 0, 1, 2. Thus, the largest finite 
group here is Z2. 

As we mentioned in the previous Section, one cannot limit oneself to subgroups of individual 
U{l)i factors or to direct products of such subgroups. In order to find all realizable groups, 
one has to write the full list of possible monomials and then calculate the symmetry group 
of all distinct pairs (for = 3) of monomials. For example, if the two monomials are vi = 
Ai(0l02)(0i03) and V2 = A2(020i)(0203)5 then the matrix X{vi + V2) defined in (!26|) has form 

^(^1 + ^'2) = ( _\ ) . (34) 

It can be diagonalized by adding the first row to the second, and then subtracting the second 
row twice from the first: 
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The diagonal for of the matrix obtained indicates that the symmetry group of the potential 
is Z3. The solution of the equation 

yields a = 27r/3 ■ A;, /3 = 0, which implies the transformation matrix of the doublets (l33ll . 

In 3HDM there are, up to complex conjugation, three bilinears and nine products of two 
bilinears transforming non-trivially under T. It is a straightforward exercise to check all 
possible pairs of monomials; in fact, by using permutations of the doublets the number of 
truly distinct cases is rather small. This procedure reveals just one additional finite group 
Z2 X Z2, which arises when at least two terms among (0|02)^, (0203)^5 (030i)^ present. 
This group is simply the group of sign flips of individual doublets. 

Thus, we arrive at the full list of subgroups of the maximal torus realizable as the symmetry 
groups of the Higgs potential in 3HDM: 

Z2, Z3, Z4, Z2 X Z2, f/(l), f/(l) X Z2, f/(l) X f/(l) . (37) 

Most of these groups were identified in [22] in the search of "simple" symmetries of 3HDM 
scalar potential. In that work a symmetry were characterized not by its group, as in our paper, 
but by providing a single unitary transformation 5* and then reconstructing the pattern in the 
Higgs potential which arises after requiring that it is S'-symmetric. In certain cases the authors 
of [22] found that the potential is symmetric under a larger group than (S), in accordance 
with the notion of realizable symmetry discussed above. 

The explicit correspondence between the seven symmetries Si, . . . , S7 of [22] and the list 
fl37|) is the following: 

^1 -> Z2 , ^2 -> U{1) realized as U{1)2 , ^3 ^ Z3 , (38) 
54 -» Z4 , ^5 ^ U{1) realized as f/(l)i , Se U{1) x Z2 , ^7 ^ U{1) x f/(l) . 

In addition to these symmetries, our list (137|) contains one more group Z2 x Z2, which was not 
found in [22] because it does not correspond to a "simple" symmetry. 

Let us also explicitly write the potentials which are symmetric under each group in (|3j 



U(l) X U(l) = T. The most general T-symmetric potential of 3HDM is given by ( 13T]) . 

U(l). This group can be realized in two non-equivalent ways, which are conjugate 
either to U{l)i or U{1)2 in ( 130|) . The distinction between the two realizations lies in the 
eigenspace decomposition: a t/(l)2-type group has a two- (complex) -dimensional space 
of eigenvectors, while eigenspaces of a t/(l)i-type group are all one-dimensional. This 
leaves more freedom in constructing a U (l)2-symmetric potential than a f/(l)i-symmetric 
one. Specifically, the general f/(l)i-invariant potential contains, in addition to f l3ip . the 
following terms: 

Al323(0l03)(403)+/^-C. (39) 

while the general ?7(l)2-invariant potential contains 



-m^3(0^03) 



-^1123(010l) + A2223(0202) + A3323(0303) 



403) 



+ A2323(0203)' + A2323(020l)(0l03) + h.C. (40) 
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It must be stressed that these potentials are written for the specific convention of groups 
U{l)i and 11(1)2 used in fl30|) . This convention reflected a specific basis used to describe 
the structure of the torus T. In other bases, the explicit terms symmetric under a 
f/(l)i or ?7(l)2-type groups will look differently. For example, the term (020i)(030i) is 
symmetric under a [/(l)i-type transformation with phases (0,a, — a). To this end, we 
note that reconstructing the symmetry group of any given potential is a separate and 
difficult task which is not addressed in the present paper (here we only classify realizable 
symmetries and give examples of their realizations). 

U(l) X Z2. Looking back to (!39l) . one might be tempted to think that the true symmetry 
group of this term is not U{1) but U{1) x Z2, because this term is also invariant under 
sign flip of 03. However, inside SU{3), a transformation with phases (0, 0, vr) is equivalent 
to (— 7r,7r, 0) which in already included in f/(l)i. 

A potential whose true symmetry group is U{1) x Z2 is given by (HOjl with all coefficients 
set to zero except for A2323- The corresponding term, (0203)^, is symmetric not only 
under the full U{1)2, but also under (— 7r,7r, 0), which generates a Z2 subgroups inside 
U{l)i. 

Z4. The potential symmetric under Z4 contains, in addition to fl3ip . the following terms: 

Al323(0l03)(403) + Ai212(0l02)' + h.C. (41) 

In accordance with the general discussion, it must contains — 1 = 2 distinct terms, as 
required for a finite symmetry group. Again, this set of terms is specific for the particular 
choice of the (f/(l)i, f/(l)2)-basis on the maximal torus. 

Z3. The potential symmetric under Z3 contains 

Al232(0l02)(0k2) + A2313(0203)(0l03) + A312I (0^0 (0^0l) + h.C. (42) 

In fact, any pair of the three terms is already sufficient to define a Z3-symmetric potential. 
Note also that different (t/(l)i, ?7(l)2)-basis on the maximal torus lead to the same Z3- 
group, because inside PSU (3) the following groups are equal: 

{{l,u,u^)) - {{co.co^,!)) - {{uj'^,l,uj)) , w = exp(27rV3) . (43) 

Z2 X Z2. This group can be realized simply as a group of independent sign flips of the 
three doublets (up to the overall sign flip). Every term in the potential must contain 
each doublet in pairs. In addition to f l3T]) . the Z2 x Z2-symmetric potential can contain 

Al212(0l02)^ + A2323(0203)^ + hM^hf + ^.C. (44) 

Z2. This group can be realized, for example, as a group of sign flips of the third doublet. 
The potential can contain any term where 03 appears in pairs. 
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4.2 The Zg x Zg-group 



The only finite abelian group that is not contained in any maximal torus in PSU{3) is Z3 x Z3. 
Although there are many such groups inside PSU{3), all of them are conjugate to each other. 
Thus, only one representative case can be considered, and we describe it with the following 
two generators 












^0 


1 






OJ 















V 





u^' ) 









/ 



exp 



27rz 
~3~ 



(45) 



A generic potential that stays invariant under this group of transformations is 



V 



—m 
+Ai 
+A3 



a. 



^03) 



2H^2) + 

(0l0l)2 + (0t02)2 + (0t03) 



+ Ao 
+ A2 



+ 



'2 + 



^^03) 



(0l02)(0l03) + {4hM<Pl) + (0^0l)(030 



h.c. 



401 P 



(46) 



with real m^, Aq, Ai, A2 and complex A3. It is also interesting to note that this group describes 
an abelian frustrated symmetry of the 3HDM potential, see details in [32] . 

The potential f H6|l is not symmetric under any continuous Higgs-family transformation, 
which can be proved, for example, using the adjoint representation of the bilinears described 
in [25] and briefly outlined in section 12.11 However it is symmetric under the exchange of any 
two doublets, e.g. 0i O 02- Together with b, such transformations generate the group S3 of 
permutations of the three doublets. Its action on elements of the Z3 generated by a is also 
naturally defined, therefore, the potential ( H6|) is automatically symmetric under the group 
(Z3 X Z3) XI Z2, which is non- abelian. 

We can conclude that the symmetry group Z3 x Z3 is not realizable for 3HDM. 



4.3 Including antiunitary symmetries in 3HDM 

In Section 13.61 we described the strategy of embedding an abelian group of unitary transfor- 
mations into a larger abelian group which includes anti-unitary transformations. Applying 
this strategy to 3HDM for all groups from the list ( 1371) . we obtain the following additional 
realizable abelian groups: 

Z; , Z2 X Z; , Z2 X Z2 X Z; , Z4 , (47) 

where the asterisk indicates that the generator of the corresponding cyclic group is an anti- 
unitary transformation. Details of the derivation are presented in Appendix [Bl Here we just 
briefly explain why groups such as Zg, Zg, U{1) x Zg do not appear in this list. When we search 
for anti-unitary transformations commuting wth the selected subgroup of T, we can indeed 
construct such groups at the price of imposing certain restrictions on the coefficients of the 
T-symmetric part (13T|) . This makes the potential symmetric under a larger group of unitary 
transformations, which is nonabelian. According to our definition, this means that the groups 
such as Zg are not realizable, although they are subgroups of larger realizable non-abelian 
symmetry groups. 
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4.4 Abelian symmetries of the 4HDM potential 

In the case of four Higgs doublets the representative maximal torus in PSU{4:) is T = f/(l)i x 
f/(l)2 X U{1)3, where 

f/(l)i = a(-l, 1, 0, 0), U{lh = (3i-2, 1,1,0), f/(l)3=^(^-^, 1 i 1^ . (48) 

The phase rotations of a generic bilinear combination of doublets under T is characterized by 
three integers p, q, r, 

i(f)l4)b) exp[i{paba + qab/3 + rabj)] (^I^b) • (49) 

These coefficients can again be represented graphically as labels of edges of three simplices 
shown in Fig. |2l 




Figure 2: Coefficients p, q, and r as labels of triangles representing bilinears in 4HDM 

With four doublets, one can construct (up to conjugation) 42 monomials transforming 
non-trivially under T. When we search for realizable finite groups, we pick up all possible 
combinations of three distinct monomials Vi, V2, W3, construct the 3x3 matrix X{vi + V2 + V3), 
check that it is non-degenerate and then diagonalize it to obtain diag((ii, d2, d^)- The symmetry 
group is then Tj^i x 1^d2 ^ ^rfg • Although this brute force algorithm can be easily implemented 
in a machine code, we can in fact apply the results of the next Section to find the list of all 
realizable finite subgroups of the maximal torus in 4HDM: 

Zfc with k = 2, ...,8; Z2 x Zfe with k = 2, 3, 4; Z2 x Z2 x Z2 . (50) 

In short, these are all finite abelian groups of order < 8. Perhaps, the most surprising of these 
groups is Z7, because it does not appear as a realizable subgroup of any single U{l)i factor 
(no two labels on Fig. [2] sum up to 7). An example of terms that possess this symmetry group 
is 

A(0l03)(0l04) + A'(40i)(404) + A"(402)(404) , (51) 

and the Z7 group is generated by the transformation with phases 

« = ^(-9, -1, 3, 7). (52) 

One can immediately check that lies in the center Z{SU{4)) and is, therefore, equivalent 
to the trivial transformation. 

Realizable continuous abelian symmetry groups as well as abelian groups containing an- 
tiunitary tansformations can also be found following the general strategy. 
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5 Abelian symmetries in general NHDM 



The algorithm described above can be used to find all abelian groups realizable as the symme- 
try groups of the Higgs potential for any A^. We do not yet have the full list of finite abelian 
groups for a generic presented in a compact form, although we put forth a conjecture 
concerning this issue, see Conjecture [6] below. However, several strong results can be proved 
about the order and possible structure of finite realizable subgroups of the maximal torus. 

Throughout this Section we will often use n := N — 1. Also, whenever we mention in this 
Section a finite abelian group we actually imply a finite realizable subgroup of the maximal 
torus. 



5.1 Upper bound on the order of finite abelian groups 

It can be expected from the general construction that for any given A^ there exists an upper 
bound on the order of finite realizable subgroups of the maximal torus in NHDM. Indeed, in 
this Section we will prove the following theorem: 

Theorem 3. The exact upper bound on the order of the realizable finite subgroup of maximal 
torus in NHDM is 

|G|<2^~^ (53) 

Before presenting the proof, let us first develop some convenient tools. First, with the 
choice of the maximal torus (12T]) . we construct n = N — 1 bilinears i = 1, . . . ,n. The 

vectors of coefficients = {pi, qi, . . . , ti) defined in fl25|) can be easily written: 

ai = (2,3,4, ...,72,1), 
02 = (1,3,4, ...,n,l), 
as = (1,2,4, ...,n,l). 



a„ = (1,2,3,. (54) 
One can use these vectors to construct the n x n matrix A: 

( a, \ 



A 



det A = 1 . (55) 



From the unit determinant we can also conclude that after diagonalization the matrix A 
becomes the identity matrix. 

Now consider a bilinear {4>\4>j) with z, j ^ 1; its vector of coefficients can be represented as 
aj_i — ai_i. More generally, for any monomial {4>\4>j){4'\4>m) with any j, /c,m, the vector of 
coefficients has the form aj_i — aj_i +am-i — ctfc-i, where oq is understood as zero. This means 
that the vector of coefficients of any monomial can be represented as a linear combination 
of a's with coefficients 0, ±1 and ±2. Since the X-matrix (l26l) is constructed from n such 
vectors, we can represent it as 

^ik = CijAjk ■ (56) 
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The square nxn matrix Cij can contain rows only of the following nine types (up to permutation 
and overall sign change): 



It follows from ( 156|) and ( 155|) that det X = det c ■ det A = det c. Therefore, order of any finite 
group is given by the module of determinant of c: 1(71 = 1 detc|. 

Let us also note two properties of the strings of type 1-9. Take any such string of length n, 
X{n) = {^1, X2, . . . , Xn), which is obtained from (1571) by an arbitrary permutation and possibly 
an overall sign flip. Then any its substring a;(„_i) = (xi, . . . , Xk-i, Xk+i, ■ ■ ■ , Xn) obtained by 
removing an arbitrary element Xk is also of type 1-9. Moreover, the element removed can be 
added at any place, and still the string x'^^_^-^ = (xi, . . . , Xm + Xk, ■ ■ ■ , Xk-i, Xk+i, ■ ■ ■ , Xn) 
remains of type 1-9. Both properties can be proved by direct inspection of all the strings. 

Now we are ready to prove Theorem [3l 

Proof. We prove by induction. Suppose that for any [n — 1) x [n — 1) square matrix -D„_i 
whose rows are strings of type 1-9, its determinant dn-i = det Dn-i is limited by |ci„_i| < 2""^. 
Take now any nxn matrix Dn constructed from the same family of strings and compute its 
determinant dn by minor expansion over the first row, with c/^^li, (ii^li, . . . , being the relevant 
minors. The procedure then depends on what type the first row is. 

1. If the first row is of type 1 or 2, then \dn\ < 2\d^^}_^\ < 2". 

2. If the first row is of type 3 or 4, then = ± < Mi-il + Mi-il < 2^^. 

3. If the first row is of type 5, then we permute the columns so that it becomes exactly 
as in fl57|) and then add the second column to the first one. Then, the first row be- 
comes (1, —1, 0, ■ ■ ■ , 0). The first minor does not change and contains rows such as 
{x2, xs, . . . , Xn), while the second minor contains rows (xi + X2, X3, . . . , Due to the 
properties discussed above, these strings are also of type 1-9, therefore the induction 
assumption applies to both minors. We therefore conclude that \dn\ = + < 2". 

4. If the first row is of type 6, then repeat follow the same procedure with only change that 
we add the second column to the third. The first row becomes (1, 1, 0, ... , 0), while 
the other rows have generic form (xi, X2, X2 + X3, X4, . . . ). The first minor contains 
rows of the form (x2, X2 + X3, X4, . . . ), which is equivalent to (x2, X3, X4, . . . ), while 
the second minor contains (xi, X2 + X3, X4, . . . ). Both rows are of type 1-9, therefore 



type 1 
type 2 
type 3 
type 4 
type 5 
type 6 
type 7 
type 8 
type 9 



(1, 0, ■■•,0), 

(2, 0, ■■•,0), 

(1,1, 0, ■■■,0), 

(1, -1, 0, ■■■,0), 

(2, -1, 0, ■■■,0), 

(1, 1, -1, 0, ■■■,0), 

(2, -2, 0, ■■■ , 0), 

(2, -1, -1, 0, ■■■,0), 

(1,1,-1,-1, 0, ■■■,0). 



(57) 



"n 




< 2". 
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5. If the first row is of type 7, we follow the procedure described for type 5 and get |(i„| = 
2|4-il <2". 

6. If the first row is of type 8, we add the second and the third columns to the first, so that 
the first row becomes (0, —1, —1, . . . , 0), while the other rows are of the form (xi + 

+ X2, X3, X4, . . . ). The second minor then contains rows {X1+X2 + X3, X3, X4, . . . ) 
which are equivalent to (xi + X2, X3, X4, . . . ), being of allowed type. The third minor 
contains {xi + X2 + X3, X2, X4, . . . ) equivalent to (xi + X3, X2, X4, . . . ), again of allowed 
type. Therefore the induction assumption applies to both minors, and we conclude that 

7. If the first row is of type 9, we add the third column to the first and the fourth column 
to the second. The first row turns into (0, 0, —1, —1, . . . , 0), while the other rows 
become of form (xi + X3, X2 + X4, X3, X4, . . . ). The third minor is built of rows (xi + 
X3, X2 + X4, X4, ... ), which are equivalent to (xi + X3, X2, X4, ... ), while the fourth 
minor is built of rows (xi + X3, X2 + X4, X3, . . . ) equivalent to (xi, X2 + X4, X3, ... ). 
Both are of the allowed type, so the induction assumption applies to both minors, and 

Therefore, \dn\ < 2" follows for any type of the first row, which completes the proof. □ 
5.2 Cyclic groups and their products 

Here we prove two propositions which show that a rather broad class of finite abelian groups 
are realizable as symmetry groups of the Higgs potential in the NHDM. First we show which 
cyclic groups are realizable and then consider direct products of cyclic groups. 

We showed in the previous subsection that in order to find the order of the finite group, 
one can use representation fl56l) and focus on the matrix Cjj instead of the matrix X. Let us 
now prove that the matrix c can also be used to find the finite symmetry group of the given 
potential. First, recall that the matrix X defines a system of linear equations on angles 

XijOj = CikAkjaj = 2'Kri . (58) 

This system remains the same upon 

• simultaneous sign change of the A;-th column in c and k-th row in A; 

• simultaneous exchange of two columns in c and two rows in A; 

• simultaneous summation of two columns in c and subtraction of two rows in A. 

• sign changes of the k-th column in A and a^; 

• exchange of columns i and j in A and exchange -H- aj; 

• summation of two columns in A and subtraction of two a's; 

• sign changes of the k-th row in c and of the integer parameters r^; 

• exchange of two rows in c and of two r's; 
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summation of two rows in c and of two parameters r. 



In short, it means that all allowed manipulations with integer matrices described earlier can 
be used for A and c. 

Then we proceeding in the following way. We first diagonalize A; the matrix c becomes 
modified by a series of allowed transformations. But diagonal A is equal to the identity 
matrix. Therefore, we are left only with c in the system of equations. We then diagonalize c 
and construct the symmetry group from its diagonal values. This means that in order to prove 
that a given group is realizable, we simply need to give an example of matrix Cij constructed 
from rows of type 1-9, fl57|) . which yields the desired diagonal values after diagonalization. 

Now we move to the two propositions. 

Proposition 4. The cyclic group Zj, is realizable for any positive integer p < 2". 
Proof. We start with the following n x n matrix Cij-. 



( 2 






C2" 



-1 



V 



2 


























2 




\ 






-1 

2 



(59) 



By straightforward manipulation with columns, starting from the last one, we arrive at 



C2" 



/ 





V 








2n-l 





2n-2 










2n-3 
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1 



(60) 



which produces the finite group Z2n. Now, let us modify by replacing the zero at the left 
lower corner by —1. Then the same transformation leads us to the group Z2n_i. If, instead, 
we replace any zero in the first column by —1, c^i = —1, then the diagonalization procedure 
leads us to the group Zp with p = 2" — 2""'^. 

Now, generally, take any integer < g < 2" and write it in the binary form. This form 
uses at most n digits. Write this binary form as a vector with n components and subtract it 
from the first column of the matrix C2»i. Then, after diagonalization, we obtain the symmetry 
group Zp with p = 2" — g. This completes the proof. □ 

Just to illustrate this construction, for n = 5 and g = 23 we write 



g=23io= IOIII2 



/ 1 \ 


1 
1 

V 1 / 



C32 - g = 



/ 1 


-1 
-1 

\ -1 



2 








-1 
2 








-1 
2 




\ 




-1 
2 / 
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Repeating the above procedure we obtain at the lowest left corner 2^ — 2^ — 2^ — 2^—2'^ = 
32 — 23 = 9, which gives the Zg group. 

Now we show that not only cyclic groups but many of their products can be realized as 
symmetry groups of some potential. 

Proposition 5. Let n = Yli=i''^i ^ partitioning of n into a sum of non-negative integers 
Hi. Then, the finite group 

G = Zp^ X X ■ ■ ■ X Zp^ (61) 

is realizable for any < j>i < 2"'. 

Proof. Let us start again from the matrix ( l59l) . For any partitioning of n one can turn this 
matrix into a block-diagonal matrix by replacing some of the — I's by zeros. The matrix is then 
represented by the smaller square blocks of size ni, 77-2, . . . , Uk- The i-th blocks has exactly 
the form of fl59l) . with n replaced by rij. Therefore, within this block one can encode any cyclic 
group Zp., with < Pi < 2"'\ Since each block can be treated independently, we can encode 
any group of the form ( 16T|) . □ 

Let us note that this proposition covers many but not all possible finite abelian groups 
with order < 2". For example, for n = 5, we can think of the group Z5 x Z5 whose order is 
smaller than 2^ = 32. However, there exists no partitioning of 5 that would lead to this group 
by applying the proposition just proved. At this moment, it remains an open problem if such 
groups are realizable. Using Theorem [3l we can formulate the following conjecture: 

Conjecture 6. Any finite abelian group with order < 2^~^ is realizable in NHDM. 

If proven true, this conjecture will give the complete classification of realizable subgroups 
of the maximal torus in NHDM. 



6 Conclusions 

Models with several Higgs doublets have become a popular framework to introduce physics 
beyond the Standard Model. When building the scalar sector of an A^-doublet model, one 
usually imposes various discrete symmetries, and many specific symmetry groups have been 
used. However, the complete lists of groups that can be implemented for a given are 
unknown for N > 2. 

In this paper we made a step towards classification of possible symmetries of the scalar 
sector of the NHDM. Namely, we studied which abelian groups can be realized as symmetry 
groups of the NHDM potential. We proved that they can be either subgroups of the maximal 
torus, or the image under the canonical homomorphism SU{N) — )■ PSU (N) of a finite nilpo- 
tent group of class 2. For the subgroups of the maximal torus, we developed an algorithmic 
strategy that gives full list of possible realizable abelian symmetries for any given N. We 
illustrated how the strategy works with two small- examples. For 3HDM we gave the full 
list of realizable abelian symmetries, including the ones with generalized CP-transformations, 
while for the 4HDM we just listed all finite realizable subgroups of the maximal torus. We 
also proved that the order of any realizable finite group in NHDM is < 2^~^ and explicitly 
described which finite groups can appear at given A^. Finally, we conjectured that any finite 
abelian group with order < 2^"^ is realizable in the NHDM. 
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It is worth stressing that the Higgs potential taken in the usual quadratic plus quartic form 
plays a major factor in restricting possible symmetry groups from all subgroups of SU{N). If 
the potential could contain higher-order terms, the list of realizable group would be larger. It 
means, in particular, that the recent impressive works on the beyond-SM applications of the 
finite subgroups of SU{3), [23], cannot be directly applied to the classification of the Higgs 
symmetries in 3HDM because no criterion was given which of the groups are realizable and 
which are not. 

An obvious direction of future research is to understand phenomenological consequences 
of the symmetries found. This includes, in particular, a study of how the symmetries of the 
potential are broken, and what is the effect of these symmetries on the physical Higgs boson 
properties. One particular result already obtained here is the possibility to accommodate 
scalar dark matter stabilized by the Zp symmetry group, for any p < 2^~^, |3]. 
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A Maximal abelian subgroups in PSU{N) 

Let V = be a unitary space of dimension N equipped with unitary form [u, v). If L is a 
subspace of V, then denotes the orthogonal complements to L in V. Now SU{N) is the 
group of all transformations of determinant 1 preserving the form, i.e., 

SU{N) = {x I detx = 1 and • X = 1}. 

We denote the group of all transformations preserving the form by GU{N). Clearly SU{N) 
is a normal subgroup of GU{N) {SU{N) < GU{N)). 

In the following lemma we collect known facts and we provide the proofs for the reader's 
convenience. 

Lemma 7. Let L be a subspace ofV andx be a linear transformation ofV. Then the following 
statements hold. 

(1) If L is x-invariant, then L'^ is -invariant. 

(2) IfxE GU{N), then x{y) = \v implies x^{v) = \v. 

(3) If X E GU{N), then there exists an orthonormal basis consisting of eigenvectors of x. 
Moreover, if Spec(x) = {Ai, . . . , A^} is the set of all eigenvalues ofx and Vx = {v E V \ 
x{v) = Xv}, then there exists the orthogonal decomposition V = Vx^ . . . Vx^. 

(4) //x, y G GU{N) commute, A is an eigenvalue ofx, then Vx is y -invariant. 
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(5) If A < GU{N) is ahelian, then A is conjugate to a subgroup of D{N)in GU{N), where 
D{N) is a group of all diagonal matrices in GU{N). In particular, if A < SU{N) is 
ahelian, then A is conjugate to a subgroup of D{N) fl SU{N) in SU{N). 

(6) Suppose H IS the stabilizer m SU{N) of L and L\ i.e., H = {GU{L)xGU{L^))nSU{N). 
Assume that L,L^ ^ 0. Let : SU{N) — )■ PSU{N) be the reduction modulo scalars. 
Then the natural projection it : H ^ GU (L) is surjective. Moreover, if T is a maximal 
torus of H, then the projection of T into GU (L) is a maximal torus of GU (L) . 

Proof. By definition of x''', for every u,v ^ V we have {x{u),v) = {u,x^{v)), whence (1). 
(2) Assume tliat x{v) = Xv. Tlien 



wlience x''"(t>) — At> = 0. 

(3) Since C is algebraically closed, there exists an eigenvalue A of x. Therefore there exists 
an eigenvector v such that x(f) = Xv. Clearly we may assume that = 1. Consider 
the line L{v) spanned by v. In view of item (2) of the lemma, L{v) is x^-invariant, hence 
item (1) of the lemma implies that L^vY is x-invariant. Therefore we obatin a decomposition 
V = L{v) © L[vy of V into an orthogonal sum of x-invariant subspaces. Induction on dim(y) 
completes the proof of this item. 

(4) For every t> G Va we have x(y(t>)) = y(x{v)) = Xy{v), i.e., y{v) G V\. Therefore V\ is a 
y-invariant subspace. 

(5) Assume that A < GU{N) is abelian. Clearly we may assume that A is nonscalar. 
Choose X e A \ Z{GU{N)). Given A G C consider Vx = {v e V \ x{v) = Xv}. Denote by 
Spec{x) the set of all eigenvalues of x. Since C is algebraically closed, we obtain that Spec(x) 
is nonempty. Moreover, since x G yl\ Z{GU{N)), it follows that Spec{x) contains at least two 
elements. In view of the proof of item (3) we obtain the orthogonal decomposition 



Choose y E A. Then, for every A G Spec(x) the subspace Vx is y-invariant, hence Vx is 
A-invariant. Induction on dim(V^) completes the proof of this item. 

(6) A matrix from {GU{L) x GU{L^)) n SU{N) has the following form: 



and l^l ■ \B\ = 1. Assume that A G GU{L) and \A\ = /i. Then = 1, i.e. /i ■ //* = 1, so if 
D = diag{^*, !,...,!), then 



= {x{v) - Xv, x{v) - Xv) = ((x - Xl){v), (x - Xl){v)) = {v, (x^ - Al)(x - Xl){v)) = 

{v, (x - Al)(xt - Xl)iv)) = ((xt - Xl){v), (xt - Al)(^)) 



V — ®\£Spec{x)Vx- 





Moreover 7r(x) = A, whence (6). 



□ 
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Consider the canonical epimorphism 



: SU{N) SU{N)/Z{SU{N)) = SU{N) ~ PSU{N) , (62) 

where Z{SU{N)) ~ Zat is the center of SU{N). Given subgroup A of PSU{N) we denote by 
A the complete preimage of A in SU{N). The following statement is evident. 

Proposition 8. Let A be an ahelian subgroup of PSU{N). Then A is either an abelian group 
or a nilpotent group of class 2 and Z{SU{N)) < Z{A). 

Proof. Since Z{SU{N)) < Z{A) for every A < SU{N), the statement is evident. □ 

If A is a maximal abelian subgroup of PSU (N) such that A is abelian, then A must be a 
maximal abelian subgroup of SU{N), and item (5) of Lemma [7] implies that A is a maximal 
torus. 

Now we find the structure of A in case, when A is nonabelian. We assume that A is a 
maximal abelian subgroup of PSU{N). 

Lemma 9. If A is nonabelian, then Z{A) = Z{SU{N)) . 

Proof Inclusion ^(^[/(iV)) < Z (A) is evident. Assume that there exists x G Z{A)\Z{SU{N)). 
Then x is nonscalar, so Spec{x) = {Ai, . . . , A^} and k > 2. In view of Lemma [7] items (3) and 

(4) , A < Csu{N)i>^) = GU{Ni) X ... X GU{Nk), where Ni = dim{VxJ. Now Lemma [7] items 

(5) and (6) implies that the projection of A into GU{Ni) is contained in a maximal torus 
of GU{Ni). Since A is a maximal abelian subgroup of PSU{N), we obtain that the pro- 
jection of A coincides with a maximal torus of GU{Ni) for every i. Hence A is a maximal 
torus of PSU{N), and so A is a maximal torus of SU{N). In particular, A is abelian, a 
contradiction. □ 

Since Z{A) = Z{SU{N)) ^ Z^r, we obtain that the exponent of Z{A) equals A^. In view 
of [331 5.2.22(i)], the exponent of A divides A^^, in particular, 7i{A) = vr(A^), where n{A) 
denotes the set of prime divisors of 1^41 and 7r(A^) denotes the set of prime divisors of A^. 

Lemma 10. A is finite. 

Proof. Let H he a. maximal normal abelian subgroup of A. Then Ga{H) = H (see [331 5.2.3]) 
and H is contained in a maximal torus of SU{N). In particular, H is generated by a finite 
number of elements. Since 7i{H) C 7r(yl) = 7r(A^) is finite and the exponent of H divides A^^, 
we obtain that H is finite. Since H = Ca{H), it follows that A/H is isomorphic to a subgroup 
of Aut{H), in particular A/H is finite. Therefore A is finite. □ 

All the information about abelian subgroups of PSU (N) obtained above was assembled in 
Theorem [T] in the main text. 

We finish by giving an example of subgroups satisfying statement (2) of Theorem [H Let 
u be the primitive A^-th root of 1 in C. Consider x = diag{l,u,u'^, . . . ,u'^~^), y = ei,2 + 
62,3 + • • • + eAr-i,Af + Cat,!, where Cij is the matrix with 1 on the (i, j)-place and on the 
remaining places. Then M = (x, y) is a nilpotent subgroup of class 2 of SU (N) for A^ odd, 
Z{N) = Z{SU{N)) and the exponent of M equals A^. 
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B Abelian symmetries in 3HDM with antiunitary trans- 
formations 



Here we explicitly describe embedding of each of the unitary abelian groups ( 137|) into a larger 
abelian group that contains antiunitary transformations. Throughout this section we will use 
the following notation: if a cyclic group TLq is generated by an antiunitary transformation, we 
indicate it by an asterisk: Z*. 



B.l Embedding U{\\ x U{\)2 

We proved in the general case that the maximal torus cannot be embedded in a larger abelian 
group. For illustration, let us repeat the argument in the specific case of 3HDM. At the first 
step of the strategy we search for a matrix h G PSU{3) that satisfies aba = b for all a & T. 
Since a is diagonal, each element {aba)ij is equal to bij up to a phase rotation ipij: 

/-2a -4/3/3 -(3/3 -a-(3/3\ 
ijjij = -13/3 2a + 2(3/3 a + 2(3/3 . (63) 

\ -a- (3/3 a + 2(3/3 2(3/3 j 

In order for an element bij to be non-zero, the corresponding phase rotation must be zero 
or multiple of 27r, modulo to 2'Kk/3 along the diagonal. It is impossible to find such b for 
arbitraty a and (3. Therefore. T = f7(l)i x f/(l)2 cannot be embedded into a larger abelian 
group. 



B.2 Embedding U{1) 

There are two distinct sorts of f/(l) groups inside T: f/(l)i-type and f/(l)2-type. Consider 
first f/(l)i. The phase rotations for the group f/(l)i {(3 = 0, arbitrary a) allow b to have 
non-zero entries 612, &21 and 633 only. We are free to choose 612 = 621 = 1, ^33 = 1 (the fact 
that det b = —1 instead of 1 is inessential; we can always multiply all the transformations by 
the overall —1 factor). Then, the transformation J' which commutes with any a G U{l)i is 



J' = 1 ■ J, af = J'a yae f/(l)i . (64) 




Next, we search for transformations x G PSU{3) such that xJ' also commutes with any 
element a G U{l)i] such transformations x form the centralizer of U{l)i. From the equation 
axa~^ = X, analyzed with the same technique of phase rotations, we find that x must be 
diagonal. By checking that the product of xiJ' and X2J' stays inside t/(l)i, we can conclude 
that X must belong to one of the ?7(l)i-type groups X^: 

x = \ e'^ I I | . (65) 

Here ^ G [0,2-7r/3) is an arbitrary but fixed parameter specifying which group we take, 
while 7 is the running angle parametrizing the elements of this group. At this stage, any 
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choice of ^ is acceptable. Thus, we embedded f/(l)i into an abehan group generated by 
(f/(l)i, J^') ~ f/(l) X Z*, where 

/ \ / \ 

J^'= \.J'=i e'^ -J. (66) 

\ e-2^« / \ e"2'« / 

Let us now see whether this group is reahzable and which must be chosen. The only 
U{l)i symmetric terms in the potential are 

A(0l03)(0|03) + \*i4^lM<p2) . (67) 

This sum of two terms is indeed invariant under J'^ provided ^ = ipx/3, where ipx is the phase 
of A. This prescription uniquely specifies which group C PSU{3). Besides, in order to 
guarantee that the T-symmetric terms (13T]) are invariant under J", we must set 

"^11 — ''^22 ' -^11 — -^22 7 -^13 = -^23 7 ^'13 = -^23 • (68) 

However upon this reduction of free parameters we observe that the resulting potential acquires 
an additional unitary symmetry: the exchange of the field two doublets 0i O 02- This 
transformation does not commute with U{l)i. Therefore, the true symmetry group of such 
potential is non-abelian, though it contains the desired abelian subgroup f/(l) x Zg. According 
to our definition, we conclude that U{1) x is not realizable in 3HDM. 

Consider now U{1)2'- a = 0, arbitrary /3. It follows from ( l63l) that no non-trivial matrix b G 
PSU{3) can satisfy aba = b. Therefore, U{1)2 cannot be embedded into a larger abelian group 
with antiunitary transformations. Note also that selecting various diagonal U{1) subgroups of 
T, for example by setting a + /3/3 = 0, will result just in another version of f/(l)i. Finally, we 
note that the unitary U{1) x Z2 emerges only when the continuous group is of the f/(l)2-type. 
Therefore, U{1) x Z2 cannot be embedded into a larger realizable abelian group as well. 

The overall result of the last two subsections is that in 3HDM all realizable continuous 
abelian groups are necessarity unitary and cannot contain antiunitary transformations. 



B.3 Embedding Z4 

The Z4 symmetry group with generator a arises as a subgroup of U{l)i. Therefore, b and 
J' can be chosen as before, (Elj), but the new conditions for x should now be checked. This 
transformation can still be represented as in fl65|) . but the angle 7 can take discrete values: 



7 = — A; or 1 = ~r ~\ — ^7 keZ. (69) 

Correspondingly, two possibilities for embedding the Z4 group arise: Z4 x Zg and Zg. 

Consider first the Z4 x Zg embedding. The Zg subgroup is generated J" in fl66|) . where, as 
before, ^ = V'a/S. As before, the potential satisfies fl68|) . includes fl67|) . but in addition it now 
contains 

A'(0^02)'+A'*(40l)^ (70) 

Note that these terms are invariant under J", and therefore, under the full Z4 x Zg. which 
is also J"-symmetric when ,^1 = ,^2- However, just as in the previous case, restrictions (l68l) 
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result in a potential which is symmetric under another Z2 group generated by the unitary 
transformation 

/ e^^ \ , 
d = e-'^ with 5 = -f , (71) 
\ 1 / ^ 

where ipx/ is the phase of A' in ( !70|) . Since d does not commute with phase rotations, therefore, 
the resulting symmetry group of the potential is non-abelian. Hence, the Z4 x symmetry 
group is not realizable. 

Consider now the Zg group generated by 

4" = I e'«+*"/4 I ■ J (72) 

e-2^« 

with the property {J"'Y = cl, the generator of the Z4. Note that upon action of J'" the term 
(0102)^ just changes its sign; therefore, such term cannot appear in the potential. But in this 
case we are left with only one type of Z4-symmetric terms. Then, according to the general 
discussion, rankXiV) = 1 and, therefore, the potential becomes symmetric under a continuous 
symmetry group. Thus, Zg is not realizable as a symmetry group of a 3HDM potential. 



B.4 Embedding Z3 

The Z3 subgroup of the maximal torus T is generated by the transformation a which rotates 
the phases of the doublets by (— 27r/3, 27r/3, 0). This group also arises from the f/(l)i-type 
group, therefore the representation for J' in 0641) and x in 0651) are still valid. Closing the 
group under product requires that 7 in fl65|) is a multiple of ir/S. We can introduce 

/ e^«~*^/3 \ 
f^= e^€+-/3 -J, (73) 

V e-2'« / 

and since {J"Y = a, this transformation generates the group Zg. 
The Zs-symmetric terms in the potential are 

Ai(0l02)(0l03) + \2{4hM<Pi) + H4<PiM<p2) + h.c. (74) 

with complex Aj. These terms become symmetric under Zg, if the following conditions are 
satisfied: 

|Ai| = IA2I and ipi + 1^2 + "ips = , (75) 

where i/ji are phases of Aj. Indeed, one can check that in this case the sum (17^ is left invariant 
under J^' with ^ = —ips/S. It is, therefore, possible to have a Zg-symmetric potential by 
requiring fl75|) and applying the usual restrictions fl68l) . 

However, just as in the previous cases, the resulting potential becomes symmetric under a 
larger unitary symmetry group, and therefore Zg is not realizable as well. The transformation 
that leaves the potential symmetric is of type 071 p with 6 = [ipi — 1^2) f'^- 
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B.5 Embedding Z2 x Z2 



The Z2 X Z2 C T symmetry can be realized as a group of simultaneous sign flips of any 
pair of doublets (or alternatively of independent sign flips of any of the three doublets). If the 
simultaneous sign flip of doublets i and j are denoted as Rij, then the group is 1, R12, R13, R23, 
with R12R13 = R23 and any {RijY = 1. This group can be embedded into Z2 x Z2 x Zj, where 
Z2 is generated by 

/ e^«i \ 
\ e-^«i-*«2 J 

At this moment, any pair ^1, ^2 can be used to construct the Z2 x Z2 x Zg group. Note also 
that this group can be also written as Z2 x Zg x Zg or Zg x Zg x Zg by redeflning the generators. 
The Z2 X Z2-symmetric potential contains the following terms 

Ai2(0l02)' + X23{4hy + A3l(40l)' + h.c. (77) 

with complex Ajj. If there phases are denoted as ipij, then the condition when ( 177|1 are 
symmetric under some Zg with an appropriate choice of .^1, ^2 is 

i^U + ^23 + ^Asi = , (78) 

or alternatively, when the product A12A23A31 is real. 

Note that in contrast to the previous cases, the full T-symmetric potential fl3ip is invariant 
under Z2 x Z2 x Zg. This guarantees that no other unitary symmetry arises in this case, and 
therefore this symmetry group is realizable. 



B.6 Embedding Z2 

The unitary abelian symmetry group Z2 can be implemented in a variety of ways, all of which 
are equivalent. One can take, for example, the Z2 group generated by sign flips of the flrst 
and second doublets, Ru- Again, the usual CP-transformation J commutes with Ru- An 
analysis similar to the previous case shows that the Z2 x Z2 group generated by Ru and J is 
realizable. 

There is also a possibility to embed the Z2 group into Z^. Consider the following transfor- 
mation 

/ \ 

4= -e*« -J, (J')2 = i?i2, (79) 
V e-2i5 J 

where ^ is an arbitrary but flxed parameter. The potential with the usual restrictions fl68l) 
and with extra terms 

Wlhf + U44>2) [(0l0i) - {4>l<p2)] + Xri^hMh) + + + h.c. (80) 

satisfying conditions [Asl = IA9I and ips + ipd = 2iljY + 7r is invariant under with 6.^ = tps + i^g. 
Adjusting A5 and Ae, one can guarantee that no other unitary symmetry appears as a result 
of (|68|) . Therefore, this Z4 group is realizable. 
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